Classical spins in topological insulators 
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Following the recent theoretical proposal and experiment on quantum spin Hall effect in 
HgTe/CdTe quantum wells, we consider a single magnetic impurity localized in the bulk of the 
system, which we treat as a classical spin. It is shown that there are always localized excited states 
in the bulk energy gap for arbitrarily strong impurity strength in inverted region, while the localized 
excited states vanish for very strong impurity strength in normal region. Similar conclusion also 
applies to three-dimensional topological insulators. This distinct difference serves as another novel 
criterion for the conventional and topological insulating phases when the time-reversal symmetry is 
broken, and can be easily experimentally observed through the STM and/or ARPES experiments. 
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PACS numbers: 75.30.Hx 73.20.At 72.25.Dc 

The quantum spin Hall (QSH) effect is a state of mat- 
ter with topological properties distinct from those of con- 
ventional insulators. PHI The first proposal of experi- 
mental realization of this effect is given in the work by 
Bernevig et al. [3] where they consider the HgTe/CdTe 
semiconductor quantum wells (QWs), and show that 
when the thickness of the QW is varied, the electronic 
states change from a normal to an inverted type at a crit- 
ical thickness d c . This transition is a topological quan- 
tum phase transition between a conventional insulating 
phase and a phase exhibiting the QSH effect with a sin- 
gle pair of helical edge states. This phase transition can 
be understood by the relativistic Dirac model in (2+1) 
dimensions, which mimic the electronic states near the T 
point. At the quantum phase transition point, d=d c , the 
mass term in the Dirac equation changes sign, leading 
to two distinct [/(l)-spin and Zi topological numbers on 
either side of the transition. Recently, the QSH phase in 
HgTe/CdTe QWs has been observed in the transport ex- 
periments, [41 which confirms Bernevig et aVs theoretical 
predictions. [3J 

Following this pioneer work, there emerge various dis- 
cussions on the novel properties of the topological in- 
sulating phase in both two- and three-dimensional (2- 
, 3-D) systems, [a-[8[ however, most of these are con- 
sidered within the framework of the preservation of the 
time-reversal symmetry (TRS), among which, we notice 
that two of them show that the novel properties of this 
topological system can also be manifested by breaking 
the TRS on the surface through the so-called topological 
magneto-electric effect Q or local charge and spin den- 
sity of states. (l(| In the meanwhile, we notice that the 
system with Mn doped impurities in the bulk of the HgTe 
QWs has been discussed in Ref.QJJ], where by breaking 
the TRS in the bulk, the quantum anomalous Hall effect 
is realized. On the other hand, it is well known that a 
single magnetic impurity in a superconductor breaks the 
TRS and induces low energy bound states in the super- 
conducting gap. [lj|, [l3| 

Motivated along this line, we discuss the presence of a 



single magnetic impurity located in the 2D bulk of the 
HgTe/CdTe QWs, which we treat as a classical spin in 
both normal and inverted regimes. Similar to the discus- 
sions in BCS superconductors by Shiba in 1968 jl3| . we 
show that in the inverted regime of the HgTe/CdTe QWs, 
there are always localized excited states (LES) in the 
bulk energy gap for arbitrarily strong impurity strength, 
while in normal regime, the LES vanish into the bulk for 
very strong impurity strength. This distinct difference 
of the response to the single magnetic impurity in bulk 
serves as another novel criteria for the conventional and 
topological insulating phases when the TRS is broken, 
and can be experimentally observed through the STM 
and/or ARPES measurements. 

The starting point of this paper is the effective four- 
band modeL3] Ho(k) = £fc<7o T o + M-k^^Tz + Ak x a^T\ + 
Aky(T T2 in HgTe/CdTe QWs around the T point in the 
basis of 1231,+), |221, +) , \E1,-), \H1,~), plus that 
of a short-range single magnetic impurity located at the 
origin. The exchange interaction in Mn doped HgTe QWs 
has been discussed in several literatures, [Tl|, [l4[ where it 
is established that the s-band and p-band electrons have 
different sp-d exchange coupling strength. To focus on 
the physical picture, we first consider the isotropic case 
where J s — J p — J, then the full Hamiltonian takes the 
form 
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Here S is the spin vector of the magnetic impurity, and 
the Pauli matrix ctj's act on spin space while Tj's act on 
the two electric subbands space, Co and r are both 2 
by 2 unit matrices. We will show later that our result is 
robust for a general form of the exchange coupling. The 
full Green's function (GF) of Hamiltonian ([T|) is obtained 
through the equation of motion formulation [131 ] as 



GW(w) = G°H<W +Gg(«)t(«)G2,M 



(2) 



where the t-matrix takes the form 

( JS^ 2 
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a F-function as 



with S = S^ + Sy + S 2 . In the above, we have introduced 



f h = jj E g 2m = dia §^ F * F - F ^ ( 4 ) 



where the diagonal elements are 



F e(h){u) = JjJ2 



m-s k + {-)M k 
D k 



(5) 



In the tight-binding model, e k = C — 2D (2 — cosfc 2 — 
cosfcy), Aik — M — 2B{2 — cosk x — cosk y ) and D k = 
(cj-£k- M k ) (w - e k + Mfc) - A 2 (sin 2 k x + sin 2 k y ), 
where M, A, B, C , D are material parameters intro- 
duced in the effective four-band model in Ref . 3] . Then 
the eigen-energies for the excited states are obtained by 
finding the poles of the GF in the bulk energy gap as 



re 

— F e(h) ( W ) = ±1 



(6) 



which is consisted of four equations. 

To characterize the momentum integration in both nor- 
mal and inverted regimes, we notice that each diagonal 
block in Ho(k) describes a quantum anomalous Hall sys- 
tem, [15| with the two forming a time-reversal conju- 
gate pair. Using the result obtained by Qi et al.,[13fc the 
topological behavior of this system is totally governed by 
two key parameters. In our case the correspondences are 
e = — Jfg + 2, which is related to the mass term in the 
(2+l)-D Dirac model, and c = — ^j- which determines 
the sign of the Chern number. Therefore in the inverted 
regime we have |e| < 2 corresponding to d > d c , other- 
wise it is topological trivial corresponding to a normal 
insulating phase. Furthermore, we set e k = without 
loss of generality to obtain a particle-hole symmetric sys- 
tem, as we know that the quadratic kinetic term has no 
contribution to the topology of this system. [15[ 

We therefore rewrite Eq.© in terms of e and c, the 
result of which with w in the bulk energy gap at several 
values of e-parameter is shown in Fig[ll where the topo- 
logical nontrival regime with |e| < 2 is given on the left 
column, while the trivial regime is plotted in the right 
column with |e| > 2. In each panel, the contribution 
from the electron subband, F e (ui), is shown in red and 
that from the hole subband, Fh(w), is shown in cyan. It is 
clear to see that the resonant condition, Eq.©, is always 
satisfied for any given impurity strength J in the inverted 
regime, and there are four LES in the bulk gap, two come 
from the electron subband and two from the hole sub- 
band. Moreover, the stronger the impurity strength is, 
the nearer the LES are to the middle of the bulk gap, 
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FIG. 1: (Color online) F-function versus cj in the bulk energy 
gap with several values of e and c 2 — I. The topological 
nontrivial regime with |e| < 2 is shown on the left column, 
while the topological trivial regime with |e| > 2 is shown on 
the right column of the figure. F e (e, w) is plotted as red lines 
while Fh(e,uj) is plotted as cyan lines. The dotted lines are 
guidelines for eyes to indicate the boundary values of the F- 
functions. 
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FIG. 2: (Color online) Energy spectrum obtained by direct 
diagonalization of Hamiltonian |T]) with eu — as a func- 
tion of impurity strength in both (a) inverted and (b) normal 
regimes. The localized states are shown in red. 
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FIG. 3: (Color online) F-function versus u) with ui in the 
bulk energy gap. The F- functions are calculated numerically 
from Eq. ([5]) with the parameters taken from Ref. J3|] for 
HgTe/CdTe QW at d = 58A and d = 70A F e (u) are shown 
as red lines while Fh(uj) are shown as cyan lines. Inset: the 
enlarged structures for the nontrivial regime. 



and no matter how strong the impurity strength is, there 
are always LES in the bulk gap which appear as peaks 
in the density of states (DOS). 

In contrast, for the normal regime, we see that the 
resonant condition is not always satisfied for any impu- 
rity strength. For weak impurity strength, there could 
exist two LES near the gap edge, however, as we in- 
crease the impurity strength, these LES merge quickly 
into the bulk and vanish finally. That is, the peaks 
in the DOS in the bulk gap are not stable against the 
strong impurity strength for this case. This distinct dif- 
ference of the response to the single magnetic impurity in 
bulk serves as another novel criteria for the conventional 
and topological insulating phases when the TRS is bro- 
ken. Note that this result is robust to the explicit form 
of the exchange interaction. For the Mn doped HgTc 
QWs [Il|, II3 |. we consider a form of exchange interac- 
tion, Jz S z a z + J* (S x <T x + S y o-y), in electron and hole 
bands separately. It turns out that the only difference 
with the isotropic model is to replace JS in Eq.© by 



(Jz' p ) 2 S 2 + (J^' P ) 2 S 2 , and in the inverted region there 
are still four persistent LES in the gap, which never merge 
into the bulk. 

To justify the above results obtained from t-matrix 
method, we directly diagonalize the Hamiltonian (JT|) on 
a square lattice in both inverted and normal regimes by 
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FIG. 4: Effective impurity concentration — ^ry as a function 
of effective impurity strength (1 - £ 2 ) 2 /(l + £ 2 ) 2 . The shaded 
area indicates the range of the impurity concentration where 
the impurity band and the continuum are separated. 



taking e = 0.5 and 3.5 for example. The obtained energy 
spectrum is plotted versus impurity strength in Fig|2] 
The four persistent LES (red lines) in the bulk energy 
gap are clearly shown for the nontrivial case where we see 
that they approach the middle of the gap as the increas- 
ing of J and do not vanish. While in the trivial regime, 
there are only two LES for small J and merge into the 
bulk for large J. This result is in perfect agreement with 
the above analysis through GF discussions. 

By using the real parameters of the HgTe/CdTe QWs 
system given in Ref. |3j, we plot the F- function at the 
quantum well width d = 58 A and d = 70 A, which 
are shown in Fig [3] We see that for d > d c where the 
QSH effect is predicted, there are always LES in the bulk 
energy gap for arbitrarily strong impurity strength; while 
for d < d c , which is a normal insulator, the LES vanish 
for strong impurity strength. 

Considering that in real materials there is always a fi- 
nite concentration of magnetic impurities, under which 
the localized excited states grow into impurity band, we 
estimate the critical concentration of magnetic impuri- 
ties -At in the topological nontrivial regime as a func- 
tion of impurity strength (1 — £ 2 ) 2 /(l + £ 2 ) 2 , at which 
the impurity band can not be distinguished from the con- 
tinuum, and the result is shown in FigSJ Here we have 
denoted £ — ^-ttN f as the effective impurity strength 
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as the effective concentration. 



However we suggest that experimentally the actual con- 
centration should be much lower than the critical values, 
so that not only its influence on the exchange coupling 
J is negligible, [IJ] but also the impurities can be con- 
sidered isolated and their coupling, such as RKKY inter- 
actions, can be ignored. Furthermore, we speculate that 
our results should be still valid even within a complete 
quantum treatment of the magnetic impurity. [l6[ 

Interestingly, the existence of LES in the bulk energy 
gap of a topological insulator is not special in two di- 



mensions, but is also true for 3D strong topological in- 
sulators (STI). As an example, we consider the strained 
HgTe which is believed to be a STI. @, [T3| The effect 
of magnetic impurities on the surface states of strained 
HgTe has been discussed by one of the authors, [HJ here 
we focus on its effect in bulk. The model Hamiltonian 
describing strained HgTe with time-reversal as well as 
inversion symmetries takes the form [{j, Il7j 



where again the results in Eqs.(j4]) and ([5]) are recovered 
[efe = automatically here since there are no kinetic 
term in the Dirac Hamiltonian (0] with A4k = Mq — 
2Mi (2 — cos k x — cos k y ) — 2M% (1 — cos k z ) and Dk = 
u 2 -[Ml +2A 2 ± (2-cosk x - cos k v ) + 2A 2 (1 -cosk z )] 
in the 3D case. The resonant conditions for LES are ob- 
tained similarly by finding the poles of the full GF in the 
bulk energy gap as 



H 3B (k) = MkT 1 + A ± k x T 5 



A ± k y T 2 



A {l k z r 3 



(7) 



where M k = M - M x (k 2 x + k 2 ) - M 2 k%, and the rep- 
resentation for Gamma matrices is chosen in such a way 
that they are invariant under the joint transformations 
of inversion and time-reversal symmetries. [18| | Two 
features are worth noticing about this model Hamilto- 
nian. Firstly, by setting A\\ = 0, Eq.([7|) recovers the 
2D HgTe/CdTe QW model. As we have discussed in 
the above, there are always LES in the inverted regime 
for arbitrarily strong impurity strength. Secondly, by 
comparing Eq.([7]) with the Kane model |19| in the ba- 
sis of \E, I) , \LH, —j) , \E, — i) , \LH, I), we observe 
that they have exactly the same form. Therefore the ma- 
trix form for Kondo-like sp-d exchange term, [ll], [lj, Il9j 
H ex (r) = — Ylm J{^ ~ Rm)S m ■ <x, in the same basis can 
be similarly extracted as 
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where A s = yN J s {S), A p = yN J p (S), [lj] and e is the 
unit vector along the direction of impurity spin. 

Following the methods developed by Fu et. al on 3D 
topological insulators with inversion symmetries, 0, 12CJ 
the STI phase characterized by odd number of Dirac 
points (k x , k y ) in the 2D surface states of Hamiltonian ([7]) 
can be analyzed through two parameters r = M2/M1 and 
n = M0/AM1. Though we wouldn't elaborate the results 
in general, some specific examples are listed below. On 
a square lattice, for r = 1.5, there is one Dirac point 
at (0,0) for r\ — 0.75; three Dirac points at (0,0), (0,7r) 
and (it, 0) for r\ = 1.25; the three Dirac points then move 
to (0,7r), (71", 0) and (n,n) for r\ — 2.25; while for n = 3 
there is only one Dirac point again at (tt,tt). For even 
larger 7*1 the system evolves out of the STI phase. 

Using the same formulation, [13| the full GF for the 
system H = J^k c \ H 3D{k)c k +^uu clH^cy is obtained 
by finding the t-matrix as 



A-A-' 



Ud (w) = 



H, 



1 - H ex F (uj) 



(9) 



A s F e (u) = ±1, 



A pV /4^3e! 



F h (u) = ±1 (10) 



By numerically plotting F e ^) & s a function of w in the 
bulk energy gap, similar behavior as shown in FigfT] is 
found respectively for STI and non-STI phases using the 
values of r and r\ illustrated above. We see that the same 
conclusion applies to 3D topological insulators. In the 
STI phase there are always LES in the bulk energy gap 
for arbitrary exchange interaction strength, while when 
out of STI phase, the LES exist only for very weak ex- 
change interaction strength. Therefore we believe that 
the existence of nonvanishing LES in the bulk energy 
gap for arbitrary impurity strength plays the role of a 
general characterization for topological insulators. 

Experimentally we suggest to detect this effect in the 
recently achieved HgTe/CdTe QWs [3, 5[ by doping a 
small concentration of Mn +2 ions in bulk. Since the 
LES evolve with the combination of exchange coupling 
strength and the magnetic moments (at the mean-field 
level), though it is hard to adjust the impurity exchange 
coupling strength, it may possible to tune the Mn mo- 
ments by a small magnetic field. [Illll4 |. When the Mn 
moments are larger than some critical value, there will 
be four peaks in the DOS spectrum in STM measure- 
ments for the QW width d > d Cl which persist for even 
larger polarization. When d < d c , the peaks in the DOS 
spectrum will be broadened and vanish as the increasing 
of the polarization. However, for 3D STI systems with 
impurities doped deep in bulk, ARPES measurements 
will be more appropriate. We suggest to use this distinct 
signal to differentiate experimentally the topological and 
the conventional insulating phases. 
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